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^ ■ After a review of the many effects of metallicity on the evolution of rotating and non-rotating 

I stars, we discuss the consequences of a high metallicity on massive star populations and on stellar 

nucleosynthesis. The most striking effect of high metallicity is to enhance the amount of mass 
lost by stellar winds. Typically at a metallicity Z = 0.001 only 9% of the total mass returned 
^ , by non-rotating massive stars is ejected by winds (91% by supernovae explosion), while at solar 

O ' metallicity this fraction may amount to more than 40%. Ifigh metallicity favors the formation 

' of Wolf-Rayet stars and of type lb supernovae. ft however disfavors the occurrence of type Ic 

supernovae. We estimate empirical yields of carbon based on the observed population of WC 
stars in the solar neighborhood, and obtain that WC stars eject between 0.2 and 0.4% of the 
mass initially locked into stars under the form of new synthesized carbon. Models give values 
well in agreement with these empirical yields. Chemical evolution models indicate that such 
carbon yields may have important impacts on the abundance of carbon at high metallicity. 
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1. General effects of metallicity on the evolution of stars 

The metallicity affects the evolution of stars mainly through its impact on the radiative 
opacities, the equation of state, the nuclear reaction rates, and the stellar mass loss rates 
(see also the discussion by Maeder 2002). Metallicity also affects the way the various 
instabilities induced by rotation occur in stars. 
O ' Recent calculations ordered by increasing metallicities are briefly presented in Table [TJ 

^ , the first column gives the reference, the second to fourth columns indicate respectively 

^ ' the range of initial masses, the values of the initial mass fraction of helium and the value 

of Z considered in the different grids of stellar models. 

In this paper we first review the various effects of metallicity on the evolution of stars, 
, then we discuss the consequences of high metallicity on the massive star populations and 

^ • their associated nucleosynthesis. 

1.1. The AF/AZ ratio 

let us recall how a change of Z , the mass fraction of heavy elements, affects the initial 
mass fraction of hydrogen, X, and of helium, Y . In the course of galactic evolution, the 
interstellar matter is progressively enriched in heavy elements and in helium. The increase 
of the abundance in helium with respect to Z is conveniently described by a t^Y j IS.Z 
law such that F = Iq + ^Y/^Z x Z (where Yq is the primordial He abundance). At 
low metallicities [Z < 0.02), Y remains practically constant, and the stellar properties 
as a function of metallicity are thus expected to be determined mainly by Z. At high 
metallicities {Z > 0.02), on the other hand, Y increases (or, alternatively, X decreases) 
significantly with Z. In those conditions, both Z and Y (and X = 1 — Y — Z) determine 
the stellar properties as a function of metallicity. 

1.2. The K-efJect 

In general, the opacities increase with increasing metallicity. This is the case for the 
bound-free (see Eq. 16.107 in Cox & Giuli 1969) and free-free (see Eq. 16.95 in Cox 
& Giuli 1969) transitions. Using the well known mass-luminosity relation L oc~ ^ ^ , 
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Table 1 . Grids of high metaUicity massive star models. 



Reference 


Initial Masses 


Y 


Z 


Claret 1997 


1 - 


40 


0.42 - 0.32 - 0.22 


0.03 


Schaerer et al. 1993 


0.8 - 


120 


0.340 


0.04 


Meynet et al. 1994 


12 - 


120 


0.320 


0.04 


Meynet & Maeder 2005^ 


20 - 


120 


0.310 


0.04 


Eldridge & Vink 2006^ 


25 - 


120 


0.340 


0.04 


Fagotto et al. 1994a 


0.6 - 


120 


0.352 


0.05 


Fagotto et al. 1994b 


0.6 


- 9 


0.475 


0.10 


Mowlavi et al. 1998a 


0.8 ■ 


- 60 


0.480 


0.10 



^ Models including the effects of rotation 

^ Models with Z-dependent mass loss rates during the WR phase. 



where fj, is the mean molecular weight, M the mass and k the opacity, one immediately 
deduces that the increase of the opacity produces a decrease of the luminosity. This can be 
seen in the left part of Fig.[T]for the 3 and 1 Mq stellar models and for metallicities inferior 
or equal to 0.04. At still higher metallicities the effect of n becomes more important than 
the opacity effect and the luminosity increases (see below). 

In contrast, for massive stars free electron scattering is the main opacity source. This 
opacity depends only on X [k^ — 0.20(1 + X)], which is about constant at Z < 0.01 
and decreases with Z at Z > 0.01. Thus at high Z the luminosity of the 20 M0 model 
shown in the left part of Fig. [T] first remains nearly constant and then increases with the 
metallicity. 

In the right part of Fig. [1] the variation of the radii on the ZAMS for stars of vari- 
ous initial masses and metallicities are shown. We see that the radii increase with the 
metallicity. From the relation L = AnR^aT^g and from the variation of L with Z, one 
can deduce that when the metallicity increases, the effective temperature decreases. For 
metallicities superior to about 0.04, the effective temperature and the luminosity both 
increase due to the ^-effect (see below). 

1.3. The ^-effect 

When the metallicity increases, ^ increases. This tends to make the star more luminous 
according to the mass-luminosity relation, provided the opacity does not increase too 
much. At very high metallicity (Z >^ 0.04), this is what happens for the 1 and 3 
models. For the 20 M0 stellar mass model, the /i-effect is also present as can be seen 
from Fig. [1] (left part). The luminosity already begins to increase as a function of Z at 
Z ~ 0.02 since the K-effect on L is negligible. 

1.4. The Enuc-ejfect 

The nuclear energy production Snuc sustains the stellar luminosity. If Enuc is arbitrarily 
increased, the central regions of the star expand, leading to a decrease in the temperatures 
and densities, and to an increase in the stellar radius. The temperature gradient therefore 
decreases leading to a decrease of both the luminosity and of the effective temperature. 

The way Enuc depends on metallicity is closely related to the mode of nuclear burning. 
When the CNO cycles are the main mode of H burning, Enuc increases with Z. In contrast, 
when the main mode of burning is the pp chain, Enuc is related to X. It is thus about 
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Figure 1. Left panels-.ZAMS locations in the HR diagram of 1, 3 and 20 Mq star models at 
different metallicities. Models ai Z — 0.1 are identified by triangles (Figure taken from Mowlavi 
et al. 1998h). Right panel: Variations of the radius as a function of initial mass, for various 
metallicities on the ZAMS. 



independent of Z at Z < 0.02, and decreases with increasing Z at Z > 0.02. These 
properties determine the behavior of Tc and pc in MS stars as a function of Z . 

Usually, the CNO cycles are the main mode of core H burning for stars more massive 
than about 1.15-1.30 Mq (depending on metallicity), while the pp chain operates at 
lower masses. At Z = 0.1, however, all stars are found to burn their hydrogen through 
the CNO cycles (at least down to M = 0.8 M© see Mowlavi et al. 1998a). 

When the CNO cycle is the main mode of burning, pc and Tc decrease with increasing 
Z (see above). This is clearly visible in the left part of Fig. [2] for the 3 and 20 M© stars at 
Z < 0.04. At Z = 0.1, however, the higher surface luminosities of the models as compared 
to those at Z — 0.04 result in concomitant higher central temperatures. 

When the pp chain is the main mode of burning, on the other hand, the location in 
the (log Pc log Tc) is not very sensitive to Z. This is well verified for the 1 M© models at 
Z < 0.1. 

The combined effects of p., k and Snuc on the stellar surface properties can be estimated 
more quantitatively with a semi-analytical approach using homology relations. This is 
developed in the appendix of Mowlavi et al. (1998b). We just recall below that a variation 
Ap of the mean molecular weight of a star of mass M affects its position in the HR 
diagram by 



AlogL.j 11 }Alog,-{ I'll }Alog.o-{ Z }Alog.o d-l) 

AiogT.,,={ }Aiog^-{ I'll }Aiog.o-{ jAiog^o (1.2) 

where eg the temperature and density independent coefficient in the relation for the 
nuclear energy production e„„c = eo T'^ (see Cox and Giuh 1969, p 692), and kq the 
opacity coefficient in the Kramers law k = kq p T~^-^. The numbers on the first line in 
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Figure 2. Left panels: Same as the left part of Fig. [T] but in the (log pc, log Tc) diagram. 
panel: Masses of convective cores relative to the stellar masses as labeled (in solar masses) next 
to the curves, as a function of metallicity. Figures taken from Mowlavi et al. (1998b). 

Eqs. 11.11 and 11.21 apply when the nuclear energy production results from the pp chain, 
while the second line applies when the CNO cycles provide the main source of nuclear 
energy. Variation of ^ is related to variation of Z by Alog/.t = j^^/iAZ. 

The mass of the convective core is mainly determined by the nuclear energy production, 
and by L. Kt Z < 0.04, it slightly decreases with increasing Z in low- and intermediate- 
mass stars, while it remains approximately constant in massive stars (see right part of 
Fig. [2). At Z > 0.04, on the other hand, it increases with Z at all stellar masseflj. We 
recall furthermore that at Z = 0.1 all stars with Al as low as 0.8 Mq possess convective 
cores. 

The MS lifetimes are summarized in the left part of Fig. [3] for several stellar masses 
as a function of metallicity. They can be understood in terms of two factors: the initial 
H abundance, which determines the quantity of available fuel, and the luminosity of the 
star, which fixes the rate at which this fuel burns. 

At Z < 0.02, tjj is mainly determined by L, being shorter at higher luminosities. The 
left part of Fig. [3] confirms, as expected from L (see left part of Fig. [T|), that increases 
with Z in low- and intermediate-mass stars, and is about independent of Z in massive 
stars. 

At Z > 0.02, on the other hand, the initial H abundance decreases sharply with 
increasing Z (the H depletion law is dictated by the adopted AY/ AZ law; for AY/ AZ = 
2.4 and Yq = 0.24, X drops from 0.69 at Z = 0.02 to 0.42 at Z = 0.1). The MS lifetimes 
are then mainly dictated by the amount of fuel available. This, combined with the higher 
luminosities at Z — 0.1, leads to MS lifetimes which are about 60% shorter at Z = 0.1 
than at Z — 0.02. This result is independent of the stellar mass for M < 40 Mq. Above 
this mass, however, the action of mass loss (see below) extends the MS lifetime, as can 
be seen from the 60 Mq curve in the left part of Fig. [31 

f It is interesting to note that we could have expected smaller convective cores in metal-rich 
massive stars where the main source of opacity is electron scattering (since this opacity is 
positively correlated with the H content). But the effect of higher luminosity in those stars 
overcomes that K-effect, and the core mass increases with increasing Z. 
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1.5. The effect of M 

When mass loss is driven by radiation, Kudritzki et al. (1989, see also Vink et al. 2001) 
showed that M in O stars is proportional to ^0-5-0-8_ As a result, the effects of mass loss 
dominate in more metal-rich stars. The stellar masses remaining at the end of the MS 
phase for various mass loss rate prescriptions are shown in the left part of Fig. [31 The 
most striking result at metallicities higher than twice solar is the rapid evaporation of 
massive stars with M > 50 Mq , and the consequent formation of WR stars during core 
H burning. 

1.6. The effects of rotation 

The effects of shellular rotation on the evolution of massive star models at high metallicity 
have been discussed in Meynet & Maeder (2005). As at lower metallicity, rotation induces 
mixing in the stellar interiors. However as was discussed in Maeder & Meynet (2001), 
rotational mixing tends to be less efficient in metal rich stars. This can be seen looking 
at the 9 M© stellar model in the left part of Fig. |4l This comes from the fact that 
the gradients of are much less steep in the higher metallicity models, so they trigger 
less efficient shear mixing. The gradients are shallower because less angular momentum 
is transported outwards by meridional currents, whose velocity scales as the inverse of 
the density in the outer layers (see the Gratton-Opick term in the expression for the 
meridional velocity in Maeder & Zahn 1998). Looking at the 40 M© stellar model in 
the left part of Fig. [4l one sees that the higher metallicity model presents the highest 
surface enrichments, in striking contrast with the behaviour of the 9 Mq model. This 
comes from the fact that in the upper part of the HR diagram the evolution of stars is 
more dominated by the mass loss rate than by rotation. And the changes occurring at 
the surface of the 40 M© are thus not only due to rotation but also to mass loss which 
is more efficient at higher Z. 

In the high mass range at high metallicity, the stars have less chance to reach the 
critical velocity during the Main Sequence phase as can be seen from the right part of 
Fig. m This is due to the fact that stellar winds remove the angular momentum at the 
surface, preventing thus the outer layers to reach the critical limit. Let us note that for 
smaller initial masses, a high metallicity may favor the approach to the critical limit. 
Indeed, in that case, the stellar winds are too weak, even at high metallicity, to remove 
a lot of mass and therefore of angular momentum, while, the meridional currents, which 
bring angular momentum from the inner regions of the star to the surface, are more rapid 
due to the lower densities achieved in the outer layers of metal rich stars. Another point 
which should be kept in mind at this point is that, at very high metallicity, the stars may 
lose large amounts of mass without losing too much angular momentum. This is due to 
the fact that in the high velocity regime, the stellar winds are polar (Maeder & Meynet 
2001). Let us stress however that this situation has little chance to be realised at high 
metallicity. Indeed the timescales for mass loss are likely shorter than the timescale for 
the torque due to wind anisotropy to affect the surface velocity. 

2. Massive star population in high metalUcity regions 

At high metallicity, as a result of the metallicity dependence of the mass loss rates, 
one expects larger fractions of Wolf-Rayet (WR) stars (Maeder et al. 1980). This can 
be seen in the left part of Fig. [5] where the WR lifetimes of rotating models for four 
metallicities are plotted as a function of the initial mass. The metallicity dependence of 
the mass loss rates is responsible for two features: 1) for a given initial mass and velocity 
the WR lifetimes are greater at higher metallicities. Typically at Z=0.040 and for M > 
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Figure 3. Left panels: Main sequence lifetimes of models of initial masses as labeled on the 
curves, as a function of their metallicity. The lifetimes are normalized for each stellar mass 
to its value a,t Z = 0.02. Right panel: Stellar masses at the end of the main sequence phase 
as a function of initial mass for different metallicitics as labeled on the curves. Dotted lines 
correspond to models computed with moderate mass loss rates, while thick lines correspond to 
models computed with twice that mass loss rates (see Mowlavi et al. 1998 for the references of 
the models used). Figures taken from Mowlavi et al. (1998b). 
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Figure 4. Left panels: Evolution during the MS phase of the N/C ratios (in number) at the 
surface of rotating stellar models as a function of the effective temperature. The differences in 
N/C ratios are given with respect to the initial values. Right panel: Evolution of the surface 
equatorial velocity as a function of time for 60 Mq stars with Wini = 300 km at different 
initial metallicities. 



60 Mq the WR lifetime is of the order of 2 Myr, while at the metallicity of the SMC the 
WR lifetimes in this mass range are between 0.4-0.8 Myr; 2) the minimum mass for a 
single star to evolve into the WR phase is lower at higher metallicity. 

Comparisons with observed populations of WR stars are shown in Meynet & Maeder 
(2005). When the variation with the metallicity of the number ratio of WR to 0-type 
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Log Mi„i/M, Log(Z/Zg) 

Figure 5. Left panels: Lifetimes of Wolf-Rayet stars from various initial masses at four different 
metallicities. All tlie models begin their evolution with Vini ~ 300 km s~^ on the ZAMS. Right 
panel: Variation of the number ratios of type Ib/Ic supernovae to type II supernovae. The crosses 
with the error bars correspond to the values deduced from observations by Prantzos & Boissier 
(2003). The dotted line is an analytical fit proposed by these authors. The continuous and 
dashed line show the predictions of rotating and non-rotating stellar models. Note that at high 
metallicity no type Ic supernovae are predicted by the models (see text). Figures taken from 
Meynet & Maeder (2002) 

stars is considered, good agreement is obtained provided models with rotation are used. 
Models without rotation predict much too small fractions of WR to 0-type stars, even at 
high metallicity. This illustrates the fact that even at high metallicity, where the effects 
of the mass loss rates are dominant, one cannot neglect the effects of rotation. 

Models well reproduce the variation with Z of the WC/WN ratio at low metallicity, 
but underestimate this ratio at high Z (see Fig. 11 in Meynet & Maeder 2005). It might 
be that the mass loss rates during the post MS WNL phase are underestimated. Also, the 
metallicity dependence of the WR stellar winds may help in resolving this disagreement 
(see Eldridge and Vink 2006). 

Current knowledge associates supernovae of type Ib/Ic with the explosion of WR stars, 
the H-rich envelope of which has been completely removed either by stellar winds and/or 
by mass transfer through Roche Lobe overflow in a close binary system. If we concen- 
trate on the case of single star models, theory predicts that the fraction of supernovae 
progenitors without H-rich envelope with respect to H-rich supernovae should be higher 
at higher metallicity. The reason is the same as the one invoked to explain the increasing 
number ratio of WR to 0-type stars with metallicity, namely the growth of the mass loss 
rates with Z. Until recently very little observational evidence has been found confirming 
this predicted behaviour. The situation began to change with the work by Prantzos & 
Boissier (2003). These authors have derived from published data the observed number 
ratios of type Ib/Ic supernovae to type II supernovae for different metallicities. The re- 
gions considered are regions of constant star formation rate. Their results are plotted in 
the right part of Fig. \5\ 

Looking at this figure, it clearly appears that rotating models give a much better fit to 
the observed data than non-rotating models. This comparison can be viewed as a check 
of the lower initial mass limit Mwne of the stars evolving into a WR phase without 
hydrogen. Let us note also at this point that at high metallicity, no type Ic supernovae 
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Figure 6. Mass of new elements ejected by stars more massive than 8 M© per mass in stars 
given by different stellar models. A Salpeter IMF has been used. The labels MM03 and MM05 
are for respectively Meynet & Maeder 2003 and Meynet & Maeder 2005. 
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Figure 7. Mass of new helium, carbon and oxygen ejected by stars more massive than 8 M0 
per mass initially in stars given by diflterent stellar models. A Salpeter IMF has been used. The 
labels MM03 and MM05 are for respectively Meynet & Maeder 2003 and Meynet & Maeder 
2005. 



are expected to occur (supernovae showing no trace of H and He) . This is a consequence 
of the high mass loss which allows the star to enter at an early stage into the WC phase 
when still a lot of helium has to be transformed into carbon and oxygen. The star thus 
keeps a high abundance of helium at its surface until the presupernova stage and explodes 
as a type lb supernova. At low metallicities, since the mass loss rates are weaker, the 
star may enter (if it does) at a later stage of the core He-burning phase, when already 
most of the helium has been transformed into carbon and oxygen (see Smith & Maeder 
1991). In that case the star may explode as a type Ic supernova. This might explain why 
long soft Gamma Ray Bursts associated to type Ic supernovae only occur in metal poor 
regions (Hirschi et al. 2006). 
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3. Massive star nucleosynthesis in high metaUicity regions 

Let us begin this section by recalling a few orders of magnitudes. When single stars 
form with a Salpeter's IMF, about 14% of the mass locked into stars consists of massive 
stars, i.e. with masses greater than 8 Mq, 25% are locked into stars of masses between 
1 and 8 Mq and 61% in stars with masses between 0.1 and 1 Mq. When all the stars 
with masses above 1 Mq have died, about 13% of the mass initially locked into stars 
is ejected by massive stars (1% remains locked into black holes or neutron stars). The 
intermediate mass stars (from 1 to 8 Mq) eject about 18.5% of the mass initially locked 
into stars (6.5% remains locked in white dwarfs). 

Figure 6 shows for various models at different metallicities the fraction of the mass in 
stars which is eventually ejected under the form of new elements. The part ejected by 
stellar winds is distinguished from that ejected at the time of the supernova explosion. 
One sees that on the whole, the fraction of the mass initially locked into stars and 
transformed into new elements by the massive stars does not depend too much on the 
model or on the metallicity. All the results are comprised between 3.5 and 4.5%. 

One notes however that more variations appear when one looks at the proportions of 
these new elements ejected by the winds and the supernova explosion. Indeed at higher 
metallicity a greater part of the new elements synthesized by stars are ejected by stellar 
winds as can be see in Fig [HI Typically the non-rotating stellar models of Schaller et al. 
(1992) predict that a stellar population at Z = 0.001 ejects in the form of new elements 
during the supernova explosion a little less than 4% of the total mass used to form the 
stars. Models at solar metallicity and higher eject about half of their new elements at 
the time of the supernova explosion and the rest through stellar winds. This may have 
a great influence on the final yields as has been shown by Maeder (1992). Indeed, mass 
loss removes matter at earlier evolutionary stages. The ejected matter has therefore been 
partially processed by nuclear burning and has a chemical composition different from the 
one it would have if the matter would have remained locked in the star. 

This effect is responsible for many specific enrichments at high metallicity: for instance, 
because of this effect, massive star models are expected to be stronger sources of "^He, 
i2q^ 22]Yg^ 26^j g^jjjj |.Q jggg important sources of ^^O at high metallicity than at 
low metallicity. The cases of helium, carbon and oxygen is shown in Fig [71 We can see 
the importance of mass loss for the carbon yields. Indeed the larger yields are obtained 
for the models computed with enhanced mass loss rates. It is interesting here to note 
that high mass loss is not a sufhcient condition for obtaining high carbon yields. High 
carbon yields are obtained only when the star enters into the WC phase at an early 
stage of the core He-burning phase. This can be seen comparing the 60 Mq stellar model 
at Z = 0.04 computed by Meynet & Maeder (2005) with the 60 Mq stellar model at 
Z = 0.02 with enhanced mass loss rate computed by Maeder (1992). These two models 
end their lifetimes with similar final masses (the model of Meynet & Maeder 2005 with 
6.7 M0 and the model of Maeder 1992 with 7.8 M©), thus the mass of new carbon ejected 
at the time of the supernova explosion is quite similar for both models and is around 
0.45 Mq. When one compares the mass of new carbon ejected by the winds, we obtain 
7.3 and 0.16 Mq for respectively the model of Maeder (1992) and the one of Meynet & 
Maeder (2005). This difference arises from the fact that the model with enhanced mass 
loss enters the WC phase at a much earlier time of the core He-burning phase, typically 
when the mass fraction of helium at the centre, is 0.43 and the actual mass of the star 
is 25.8 Mq, while the model of Meynet & Maeder (2005) enters the WC stage when 
is 0.24 and the actual mass 8.6 Mq. 

The example above shows that the entry at an early stage into the WC phase is more 
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favored by strong mass loss than by rotation. This comes from the fact that rotation 
favors an early entry into the WN phase, while the star has still an important H-rich 
envelope. It takes time for the whole H-rich envelope to be removed and when it is done, 
the star is already well advanced into the core He-burning phase. Of course this conclusion 
is quite dependent on the magnitudes of the mass loss rates. For instance higher mass 
loss rates during (a part of) the WNL phase would favor an early entrance into the WC 
phase. This would give a better agreement with the number ratio of WC to WN observed 
and would lead to higher carbon yields. 

Let us conclude this paper by estimating an empirical yield in carbon from the Wolf- 
Rayet stars in the solar neighborhood. From the catalogue by van der Hucht (2001) one 
obtains that the number of WC stars at a distance less than 3 kpc around the Sun is 
44. The mass loss rate during the WC phase is estimated to be 10^^'^ M© per year, the 
mass fraction of carbon observed in the WC stellar wind is around 0.35 (see Table 2 in 
the review by Crowther 2006). If we consider a star formation rate of 2-4 Mq per square 
parsec and per Gyr, one obtains that the mass of (new) carbon ejected by WC wind per 
mass used to form stars is between 0.25 and 0.5%. These empirical yields are well in lines 
with the range of values given by the models which are between 0.2-0.6% (see Fig. [3 the 
greatest value corresponding to the case of Maeder 1992). Let us note that incorporating 
the yields of Maeder (1992) into chemical evolution models has an important impact (see 
e.g. Prantzos et al. 1994). If the upper value of the empirical yields is the correct one 
then this indicates that WC stars are very important sources of carbon in metal rich 
regions. 
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